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ABSTRACT 

The recent detection of Sagittarius A* at A = 1.3 mm on a baseline from Hawaii to Arizona demonstrates 
that millimeter wavelength very long baseline interferometry (VLBI) can now spatially resolve emission from 
the innermost accretion flow of the Galactic center region. Here, we investigate the ability of future millimeter 
VLBI arrays to constrain the spin and inclination of the putative black hole and the orientation of the accretion 
disk major axis within the context of radiatively inefficient accretion flow (RIAF) models. We examine the 
range of baseline visibility and closure amplitudes predicted by RIAF models to identify critical telescopes for 
determining the spin, inclination, and disk orientation of the Sgr A* black hole and accretion disk system. We 
find that baseline lengths near 3 GA have the greatest power to distinguish amongst RIAF model parameters, 
and that it will be important to include new telescopes that will form north-south baselines with a range of 
lengths. If a RIAF model describes the emission from Sgr A*, it is likely that the orientation of the accretion 
disk can be determined with the addition of a Chilean telescope to the array. Some likely disk orientations 
predict detectable fluxes on baselines between the continental United States and even a single 10-12 m dish in 
Chile. The extra information provided from closure amplitudes by a four-antenna array enhances the ability of 
VLBI to discriminate amongst model parameters. 

Subject headings: black hole physics — accretion, accretion disks — submillimeter — Galaxy: center — 
techniques: interferometric 



1. INTRODUCTION 

The Galactic center radio source Sagittarius A* is believed 
to be associated with a black hole with a ma ss of approxi- 
mately 4 x 10 6 MfT, at a distance of about 8 kpc (Schod el et all 
12003b iGhez et al.ll2008h . The event horizon of Sgr A* has 
the largest apparent angular size of all known black holes as 
viewed from Earth. High angular resolution is critical to un- 
derstanding Sgr A* because the size scales are so small: for 
instance, the apparent lensed horizon size is 55 /ias. The quest 
for angular resolution has driven observers toward very long 
baseline interfe rometry (VLBI) at mil limeter wavelengths. 
Most recently, iDoeleman et all d2008al) detected Sgr A* at 
230 GHz (A = 1.3 mm) on baselines between the James Clerk 
Maxwell Telescope (JCMT) in Hawaii and the Submillime- 
ter Telescope Observatory (SMTO) in Arizona, as well as 
between the SMTO and a Combined Array for Research in 
Millimeter-wave Astronomy (CARMA) telescope in Califor- 
nia. The former is the longest baseline (3.5 GA, fringe spacing 
~ 60 /ias) on which Sgr A* has ever been detected. 

In the absence of a true ab initio simulation of Sgr A*'s ac- 
cretion flow, a variety of physically simplified models have 
been proposed. Constraints from the observed spectrum, po- 
larization, and variability have not proved sufficient to conclu- 
sively establish the nature of the emission region. This is ev- 
idenced by the continuing vigorous debate over the morphol- 
ogy of this region, e.g., primar ily from the inner most regions 
of an inefficient accretion flow (Yuan et alj2003l) . the acceler- 
ation region of a nasce nt iet dFalcke et al . 2000), or something 
qualitatively different dlgumenshchev et alJ.12003!) . However, 
millimeter VLBI's ability to resolve the emitting region di- 
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rectly promises to discriminate amongst these. To do so will 
require a detailed understanding of the role that key physical 
parameters (such as black hole spin, disk/jet orientation, etc.) 
play in shaping the millimeter images of Sgr A*. Conversely, 
if the context of the emission can be conclusively established, 
millimeter VLBI has great potential to extract these parame- 
ters. 

In a rela ted work.lBroderick et all d2008f) examine the abil- 
ity of the IDoeleman et al.l d2008al) VLBI detections to dis- 
criminate amongst various radiatively inefficient accretion 
flow (RIAF) models parameterized by black hole spin, view- 
ing angle, and disk orientation. A key finding is that the de- 
tected flux on the JCMT-SMTO baseline argues against ori- 
entations in which the assumed accretion disk is close to face 
on. However, the RIAF parameters are strongly coupled, and 
there are not yet enough millimeter VLBI measurements to 
make strong statements about the black hole spin, for in- 
stance. In this work, we explore the question of which future 
millimeter VLBI observations are likely to have the great- 
est impact for distinguishing between RIAF models and are 
therefore most likely to determine the physical parameters of 
the black-hole/accretion-disk system. 

2. DISK MODELS AND TELESCOPES 

Here we consi der an ensemble o f RIAF models similar to 
that described in Yuan et all d2003l) and described in detail in 
IDoeleman etafl d2008bl) and iBroderick et al.l d2008h . These 
are generally characterized by the inefficient energy transfer 
between the electrons and the ions. In particular, we model 
the millimeter flux from Sgr A* as synchrotron emission due 
to populations of thermal (though with considerably lower 
temperatures than the ions) and nonthermal electrons in a ge- 
ometrically thick, quiescent disk containing a toroidal mag- 
netic field. The thermal and nonthermal electron densities and 
thermal electron temperature are solved for using measured 
radio, submillimeter, and near-infrared fluxes as constraints. 
Model images as viewed from the Earth are produced using 
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fully general relativistic radiative transfer. 

Each model has a total flux density of 2.4 Jy at A = 
1.3 mm, in line with the integrated flux density measured by 
iDoeleman et al.l d2008al) . These models are parameterized by 
black hole spin (a), inclination (9 = 0° for a face-on disk), and 
orientation of the disk major axis on the plane of the sky (£) 
measured counterclockwise from east (equivalently the posi- 
tion angle of the spin axis). (Henceforth we shall use the word 
"models" to refer to RIAF models with these three free pa- 
rameters.) Models are produced at spins from a = 0.0 to 0.9 
in increments of 0.1, as well as at 0.99 and 0.998. The en- 
semble includes models at 9 = 1° and then at 9 = 10° to 90° 
in 10° increme nts. This particular ens emble is described in 
greater detail in Broderic k et al.l d2008). We consider models 
with 0° < £ < 180° with an increment of A£ = 5° (in £□) or 
30° (in § jj3.2| and |3.3l l. Visibility amplitudes are symmetric 
under 180° rotation for any real intensity distribution. 

We consider seven potential stations for millimeter VLBI: 
Hawaii, consisting of one or more of the JCMT, Submil- 
limeter Array (SMA), and Caltech Submillimeter Observa- 
tory (CSO) possibly phased together into a single aperture; 
the SMTO on Mount Graham, Arizona; CARMA telescopes 
in California, either individually or phased together; the Large 
Millimeter Telescope (LMT) on Sierra Negra, Mexico; a 
Chilean station consisting of either the Atacama Submillime- 
ter Telescope Experiment (ASTE), Atacama Pathfinder Ar- 
ray (APEX), or a phased array of Atacama Large Millimeter 
Array (ALMA) dishes; the Institut de radioastronomie mil- 
limetrique (IRAM) 30-m dish on Pico Veleta (PV), Spain; and 
the IRAM Plateau de Bure (PdB) Interferometer, phased to- 
gether as a single aperture. Assumed telescope capab ilities 
and sensitivities are detailed in Doele man et 
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3. BLACK HOLE PARAMETER ESTIMATION 

3.1. Favored Regions of the (u,v) Plane 

For the purpose of distinguishing different models, it is 
most desirable to obtain observations of Sgr A* on baselines 
at points at which the visibilities are most disparate. We com- 
pute model visibilities 



V(u,v)= / / dxdyl(x,y)e 



-2m(xu+yv)/\ 



where I(x,y) is the model intensity distribution with x and 
y are aligned east and north, respectively. We characterize 
the variation in model visibilities via the standard deviation 
a(u, v) amongst all models in the ensemble. While we initially 
ignore the effects of interstellar scattering, it can be included 
by multiplying the visibilities by a unit-normalized elliptical 
Gaussian centered at the origin (half width at half maximum 
of 7.0 x 3.8 GA, p osition angle 170° east of north, based on 
iBower et al.ll2006h . 

Figure Q] shows the scatter in predicted model visibilities, 
a(u,v). The peak standard deviation occurs at a baseline 
length of approximately 3 GA at A = 1.3 mm (~ 4000 km) 
aligned with the accretion disk major axis. Since the disk 
orientation (£) is not known a priori, the optimal baseline ori- 
entation is unknown. 

Not all models are equally likely, however. The 
IDoeleman etafl (l2008al) detections on the JCMT-SMTO and 
SMTO-CARMA baselines as well as total flux density mea- 
surement with the CARMA interferometer place constraints 
on RIA F models. Following the approach of [Broderick et al.l 
(2008), we com pute the probability p( a, 9,£) for every model 
based on the IDoeleman et al.l (l2008al) detections, including 




FIG. 1. — Map of standard deviation of predicted RIAF model visibilities 
(cf(u,v)) at A = 1.3 mm holding £ = 0° constant. Units are in Jy. The largest 
spread in predicted model quantities occurs when the projected baseline is 
aligned with the disk major axis. Rotating the image on the plane of the sky 
produces a rotation of the visibility amplitudes in the (u, v) plane. 




FIG. 2. — Plot of the standard deviation <j p (u, v) restricted to the most proba- 
ble models (see j|3.H on the inner part of the (u, v) plane at A = 1 .3 mm. Units 
are in Jy. Potential (u. v) tracks due to Earth rotation are superposed, with let- 
ters indicating antennas. A: APEX/ASTE/ALMA, B: PdB, C: CARMA, H: 
Hawaii (incl. JCMT), L: LMT, S: SMTO. V: PV 

a sin 9 prior on the inclination. Restricting the ensemble of 
models to those with p > 0.01 p mdx (or approximately one- 
quarter of the total set of models), we can obtain the scatter 
in the predicted correlated flux densities of the most probable 
models, a p (u,v) (Fig. 13. We find that the largest scatter in 
models still occurs at baseline lengths near 3 GA, but at orien- 
tations substantially different from that of the JCMT-SMTO 
detection (labelled "HS" in Fig. [2]), which has already placed 
a stringent constraint on models in its region of (u, v) space. 

3.2. Implications for Specific Baselines 

Figure[3]shows predicted correlated flux densities on poten- 
tial millimeter VLBI baselines for the probable RIAF mod- 
els. Each pair of telescopes produces visibility measurements 
along a range of projected baseline lengths due to the Earth's 
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FIG. 3. — Predicted model fluxes on potential baselines. Each panel shows the predicted track of correlated flux density with V u 2 + v 2 (changing due to Earth 
rotation) for one baseline for all models (A£ = 30°) for which p > 0.01 p max - Lines show the maximum and minimum predicted correlated flux density from the 
set of probable models with a given disk orientation, shown in color. The visibility amplitudes on many baselines are highly sensitive to £. Baselines between 
existing VLBI stations in the continental US and new stations in Chile and Mexico (LMT) are optimal for further constraining RIAF model parameter space. 



rotation. These correlated flux densities have been corrected 
for interstellar scattering. 

Detections have already been obtained on the Hawaii- 
SMTO and SMTO-CARMA baselines, although Sgr A* 
was not detected on the Hawaii-CARMA baseline 
dDoeleman et al.1 l2008al) . It is likely that the next mil- 
limeter VLBI observations of Sgr A* will occur on an array 
containing these same three stations, although with the inclu- 
sion of a phased-array processor at Hawaii to sum the signals 
from the JCMT, CSO, and SMA. The increased sensitivity 
should allow Sgr A* to be detected on the Hawaii-SMTO 
baseline with a higher signal-to-noise ratio, tightening the 
constraints on RIAF model parameters especially if Sgr A* 
can be detected on scans when it is low to the horizon from 
Arizona, where the range of predicted model flux densities is 
larger. The Hawaii-CARMA baseline will be more sensitive 
than the Hawaii-SMTO baseline once the planned array 
phasing capability is added to CARMA, although the range 
in predicted model flux densities is somewhat smaller on 
Hawaii-CARMA. Because the predicted flux densities on the 
Hawaii-CARMA and Hawaii-SMTO baselines are similar 
due to their location in approximately the same location of 
the (u,v) plane, strongly constraining key black hole and 
accretion disk parameters will necessitate obtaining flux 
density measurements on other baselines as well. However, 
measurements on the Hawaii-CARMA baseline may be 
important to establish the validity of RIAF models of Sgr A* 
in general and rule out simpler, less physically-motivated 
models, such as a Gaussian or a ring. 

The optimal baselines for further constraining model pa- 
rameter space are SMTO-Chile, CARMA-Chile, and LMT- 
Chile, which are in the preferred zone of baseline lengths 
and at a position angle nearly orthogonal to that of Hawaii- 
SMTO. Many of the models predict flux densities of well 
over 500 mJy (some even above 1 Jy) on the SMTO-Chile 
and CARMA-Chile baselines. These flux densities should be 
detectable even if the Chilean station consists of a single 10 
or 12 m telescope, such as APEX, ASTE, or a single ALMA 
dish. 

Of the three parameters which we consider in this work, 
millimeter VLBI will most easily be able to constrain the ori- 
entation of a RIAF disk. Disk emission comes predominantly 



from the Doppler-brightened approaching side and is elon- 
gated perpendicular to the projected disk major axis when in- 
clined relative to the line of sight. Because the locations of 
existing and potential millimeter telescopes define preferred 
directions on the sky for measuring this emission, baseline 
correlated flux densities depend strongly on disk orientation. 
The scatter within the predicted correlated flux density for 
each baseline is much smaller when restricted to one value of 
£ than when the entire model parameter space is considered 
as a whole (Fig.|3). 

Constraining the disk orientation of Sgr A* may be a nec- 
essary first step toward extracting the spin and inclination of 
the system. For a given orientation, increasing the black hole 
spin or the inclination of the system causes the emitting re- 
gion to appear smaller, raising the correlated flux density on 
longer baselines (Fig. |4). Indeed, if the disk orie ntation can 
be determined, the existing Doele man et al.l (l2008al) data may 
already place strong c onstraints on the spin and inclination 
BrodericketaD (12008|). 



3.3. Potential Complications 

The flux density of 2.4 Jy at 230 GHz detected by 
iDoeleman et alJ (|2008al) is low compare d to measurements 
with the SMA dMarrone et al.l2007ll2008l) and high compared 



to me asurements with BIMA (IZhao et al.ll2()03HBower et al.l 
120051) and a PV-PdB VLBI measurement at 215 GHz 
dKrichbaum et al.ll 1 998b . Sgr A* exhibits variability, although 
the consistent fluxes m easured on two consecutive days by 
IDoeleman et al.l (l2008al) suggest that Sgr A* was observed in 
its quiescent state. Since the variability mechanism is not un- 
derstood, it remains unclear how variability in the total flux 
density of Sgr A* at millimeter wavelengths will affect the 
correlated flux densities measured on diverse VLBI baselines. 
Multiple epochs of observation will be essential to explor- 
ing the link between total flux variability and spatial structure 
around Sgr A*. It will also be important to obtain contem- 
poraneous measurements of fluxes on a variety of baselines 
as well as a simultaneous "zero-spacing" flux density, prefer- 
ably from a connected-element interferometer. 

Variable and elevation-dependent antenna gains may limit 
the accuracy to which telescope flux scales can be derived, in- 
troducing systematic errors into calculations of correlated flux 
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FIG. 4. — Left three panels: Correlated baseline flux density plots on baselines bet ween Hawaii, SMTP, and CARMA for a likely model (a = 0.6,6 = 70°, £ = 
30°), varying one parameter at a time. Data from the second day of observations by Doeleman et al. (2008a) are superposed. Increasing the black hole spin and 
inclination predict greater flux densities on longer baselines. Rotation of the disk emission in the plane of the sky can cause the predicted flux density to increase 
or decrease depending on the baseline orientation. Right panel: The same for six different likely models showing the widely divergent predicted flux densities on 
the CARMA-Chile baseline. 



densities. These systematic errors may be difficult to charac- 
terize accurately and may limit the ability of baseline-based 
quantities to distinguish between different models. Closure 
amplitudes are robust against station-based gain variations 
and will provide model constraints with significantly reduced 
systematic errors. If observations are taken with at least four 
telescopes, closure amplitudes can be constructed from ratios 
of baseline visibilities (A a b c d = \V a b\ \Vcd\ \Vac\~ 1 |Vm| -1 , where 
subscripted letters identify telescopes). 

The models predict a wide range of closure amplitudes on 
quadrangles of four telescopes, as shown in Figure [5] While 
even a single measurement of a closure amplitude will be use- 
ful, multiple measurements of closure amplitudes on a quad- 
rangle with time will have even greater power in further con- 
straining model parameter space, since different models can 
predict strongly differing closure amplitudes over a night of 
observations. As is the case for baseline visibilities, closure 
amplitudes are most highly sensitive to variations in £ among 
the most probable models. From a closure amplitude perspec- 
tive, either a Chilean telescope or the LMT would make an 
excellent fourth station in an observing array, although the 
low predicted correlated flux densities on the Hawaii-Chile 
baseline may make measurement of a second closure ampli- 
tude on an array consisting only of Hawaii, SMTO, CARMA, 
and Chile difficult without the use of phased- ALMA as the 
Chilean station. 

An alternative observational strategy in the absence of a 
fourth telescope site would be to use two different telescopes 
at the same location, for instance at the CARMA site or in 
Hawaii. A single nontrivial closure amplitude can be formed 
from the resulting visibilities. The closure amplitude result- 
ing from two CARMA antennas, the SMTO, and Hawaii has 
less power to differentiate among models than other four- 
antenna arrays (Fig. 0, but multiple detections over a sin- 
gle scan may still be useful. The redundant information pro- 
vided by the double baselines between the CARMA site and 
another antenna will also aid in the detection of fringes and 
reduce the systematic uncertainty in correlated flux densities 
on the CARMA baselines. Given th e nondetection of Sgr A * 
on the JCMT-CARMA baseline by iDoeleman et al.l (l2008al) . 
it is likely that this strategy would require the use of a phased- 
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FIG. 5. — Predicted model closure amplitudes for all probable RIAF mod- 
els over the time range for which Sgr A* is above 5° elevation at at least four 
of the five Western hemisphere stations. Colors are as in Figure fS] Small 
jumps in the Hawaii-SMTO-Chile-CARMA panel are artifacts in precision 
when the Hawaii-Chile baseline passes through a null. A single closure am- 
plitude can also be formed by using two telescopes at the same site, such as 
CARMA. Closure amplitudes are highly sensitive to £ and are robust against 
flux calibration errors. 

array processor at Hawaii. 

4. CONCLUSIONS 

The planned evolution of millimeter VLBI capability over 
the next few years will place strong constraints on RIAF mod- 
els of Sgr A* emission. It is important to devise an observ- 
ing strategy to test the RIAF hypothesis and extract informa- 
tion on the parameters of the putative accretion disk system of 
Sgr A*, including the disk orientation, inclination, and black 
hole spin. By identifying differences between models that are 
consistent with present VLBI observations, we identify four 
key points to guide future millimeter VLBI observations: 

1 . There is a preferred zone of baseline lengths near 3 GA 
where there is enough angular resolution to distinguish be- 
tween many of the models but not so much resolution that 
a large amount of flux is lost. The JCMT-SMTO detections 
are in this zone. The construction of a phased-array proces- 
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sor for the Hawaiian millimeter telescopes will allow for even 
tighter constraints to be placed on RIAF model parameters on 
the Hawaii-SMTO baseline as well as increase the chance of 
detecting Sgr A* on other baselines to Hawaii. 

2. There is a strong correlation between fluxes on the 
Hawaii-SMTO and Hawaii-CARMA baselines, limiting the 
ability of the latter baseline to place extra constraints on 
RIAF parameters unless a phased-array processor is installed 
at CARMA, making Hawaii-CARMA a more sensitive base- 
line than Hawaii-SMTO. Nevertheless, the Hawaii-CARMA 
measurement will be important both for ruling out simpler 
models (e.g., a Gaussian or annulus) of the emission from 
Sgr A* that the present VLBI measurements cannot and for 
testing the RIAF hypothesis. 

3. VLBI arrays that include either the LMT or a Chilean 
telescope provide the best constraints on RIAF models. Pre- 
dicted flux densities on the SMTO-LMT and CARMA-LMT 
baselines are well in excess of 1 Jy. Thus, the LMT can sig- 
nificantly enhance Sgr A* VLBI even with a surface not ad- 
justed for maximum accuracy. Likewise, since many models 
predict flux densities on the SMTO-Chile and CARMA-Chile 
baselines well in excess of 0.5 Jy, inclusion of even a single 



10 m or 12 m Chilean telescope will result in important model 
constraints on Sgr A*. 

4. Closure amplitudes will be critical, as amplitude calibra- 
tion errors impose systematic errors on measured baseline flux 
densities that reduce their ability to constrain RIAF model 
parameters. At least four telescopes should be used in an 
observing array in order to obtain closure amplitudes. If no 
fourth telescope site is available at the time of reobservation, 
a second antenna at the CARMA site should be included in a 
A = 1 .3 mm observing array. 

We note the important caveat that our analysis assumes 
that a RIAF accurately and completely describes the emis- 
sion from Sgr A* at 230 GHz. This assumption may be in- 
correct and is likely at best an approximation to the actual, 
complex emitting region in Sgr A*. For these reasons, it will 
be important to eventually obtain data on as many different 
millimeter-wavelength VLBI baselines as is feasible. 



The high-frequency VLBI program at Haystack Observa- 
tory is funded through a grant from the National Science 
Foundation. 
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